Search for: All records

Abstract The marine atmospheric boundary layer (ABL) and oceanic boundary layer (OBL) are a two-way coupled system. At the ocean surface, the ABL and OBL share surface fluxes of momentum and buoyancy that incorporate variations in sea surface temperature (SST) and currents. To investigate the interactions, a coupled ABL–OBL large-eddy simulation (LES) code is developed and exercised over a range of atmospheric stability. At each time step, the coupling algorithm passes oceanic currents and SST to the atmospheric LES, which in turn computes surface momentum, temperature, and humidity fluxes driving the oceanic LES. Equations for each medium are time advanced using the same time step but utilize different grid resolutions: the horizontal grid resolution in the ocean is approximately four times finer, e.g., (Δx_o, Δx_a) = (1.22, 4.88) m. Interpolation and anterpolation (its adjoint) routines connect the atmosphere and ocean surface layers. In the simplest setup of a statistically horizontally homogeneous flow, the largest scale ABL turbulent shear-convective rolls leave an imprint on the OBL currents in the upper layers. This result is shown by comparing simulations that use coupling rules that are applied either instantaneously at everyx–ygrid point or averaged across anx–yplane. The spanwise scale of the ABL turbulence is ∼1000 m, while the depth of the OBL is ∼20 m. In these homogeneous, fully coupled cases, the large-scale spatially intermittent turbulent structures in the ABL modulate SST, currents, and the connecting momentum and buoyancy fluxes, but the mean profiles in each medium are only slightly different. 

Abstract Langmuir turbulence, a dominant process in the ocean surface boundary layer, drives substantial vertical mixing that influences temperature, salinity, mixed layer depth, and biogeochemical tracer distributions. While direct resolution of Langmuir turbulence in ocean and climate models remains computationally prohibitive, its effects are commonly parameterized, frequently within established turbulent mixing frameworks like the K‐profile parameterization (KPP). This study utilizes a modified KPP that determines boundary layer depth through an integral criterion, diverging from the conventional KPP's dependence on the bulk Richardson number. The modified KPP demonstrates markedly lower sensitivity to model vertical resolution than its conventional counterpart. Building upon this modified KPP framework, we introduce an innovative parameterization scheme for Langmuir mixing effects. We evaluate the performance of this new scheme against existing approaches using a one‐dimensional (1D) column model across four different scenarios, incorporating validation against both large eddy simulation (LES) results and field measurements. Our analysis reveals that the new Langmuir mixing scheme, explicitly designed for the modified KPP framework, performs competitively while maintaining reduced sensitivity to vertical resolution. 

Abstract Realistic computational simulations in different oceanic basins reveal prevalent prograde mean flows (in the direction of topographic Rossby wave propagation along isobaths; aka topostrophy) on topographic slopes in the deep ocean, consistent with the barotropic theory of eddy-driven mean flows. Attention is focused on the western Mediterranean Sea with strong currents and steep topography. These prograde mean currents induce an opposing bottom drag stress and thus a turbulent boundary layer mean flow in the downhill direction, evidenced by a near-bottom negative mean vertical velocity. The slope-normal profile of diapycnal buoyancy mixing results in downslope mean advection near the bottom (a tendency to locally increase the mean buoyancy) and upslope buoyancy mixing (a tendency to decrease buoyancy) with associated buoyancy fluxes across the mean isopycnal surfaces (diapycnal downwelling). In the upper part of the boundary layer and nearby interior, the diapycnal turbulent buoyancy flux divergence reverses sign (diapycnal upwelling), with upward Eulerian mean buoyancy advection across isopycnal surfaces. These near-slope tendencies abate with further distance from the boundary. An along-isobath mean momentum balance shows an advective acceleration and a bottom-drag retardation of the prograde flow. The eddy buoyancy advection is significant near the slope, and the associated eddy potential energy conversion is negative, consistent with mean vertical shear flow generation for the eddies. This cross-isobath flow structure differs from previous proposals, and a new one-dimensional model is constructed for a topostrophic, stratified, slope bottom boundary layer. The broader issue of the return pathways of the global thermohaline circulation remains open, but the abyssal slope region is likely to play a dominant role. 

Estuaries, as connectors between land and ocean, have complex interactions of river and tidal flows that affect the transport of buoyant materials like floating plastics, oil spills, organic matter, and larvae. This study investigates surface-trapped buoyant particle transport in estuaries by using idealized and realistic numerical simulations along with a theoretical model. While river discharge and estuarine exchange flow are usually expected to export buoyant particles to the ocean over subtidal timescales, this study reveals a ubiquitous physical transport mechanism that causes retention of buoyant particles in estuaries. Tidally varying surface convergence fronts affect the aggregation of buoyant particles, and the coupling between particle aggregation and oscillatory tidal currents leads to landward transport at subtidal timescales. Landward transport and retention of buoyant particles is greater in small estuaries, while large estuaries tend to export buoyant particles to the ocean. A dimensionless width parameter incorporating the tidal radian frequency and lateral velocity distinguishes small and large estuaries at a transitional value of around 1. Additionally, higher river flow tends to shift estuaries toward seaward transport and export of buoyant particles. These findings provide insights into understanding the distribution of buoyant materials in estuaries and predicting their fate in the land–sea exchange processes. 

Abstract Recent years have seen a surge in interest for leveraging neural networks to parameterize small-scale or fast processes in climate and turbulence models. In this short paper, we point out two fundamental issues in this endeavor. The first concerns the difficulties neural networks may experience in capturing rare events due to limitations in how data is sampled. The second arises from the inherent multiscale nature of these systems. They combine high-frequency components (like inertia-gravity waves) with slower, evolving processes (geostrophic motion). This multiscale nature creates a significant hurdle for neural network closures. To illustrate these challenges, we focus on the atmospheric 1980 Lorenz model, a simplified version of the Primitive Equations that drive climate models. This model serves as a compelling example because it captures the essence of these difficulties. 

Abstract Release of iron (Fe) from continental shelves is a major source of this limiting nutrient for phytoplankton in the open ocean, including productive Eastern Boundary Upwelling Systems. The mechanisms governing the transport and fate of Fe along continental margins remain poorly understood, reflecting interaction of physical and biogeochemical processes that are crudely represented by global ocean biogeochemical models. Here, we use a submesoscale‐permitting physical‐biogeochemical model to investigate processes governing the delivery of shelf‐derived Fe to the open ocean along the northern U.S. West Coast. We find that a significant fraction (∼20%) of the Fe released by sediments on the shelf is transported offshore, fertilizing the broader Northeast Pacific Ocean. This transport is governed by two main pathways that reflect interaction between the wind‐driven ocean circulation and Fe release by low‐oxygen sediments: the first in the surface boundary layer during upwelling events; the second in the bottom boundary layer, associated with pervasive interactions of the poleward California Undercurrent with bottom topography. In the water column interior, transient and standing eddies strengthen offshore transport, counteracting the onshore pull of the mean upwelling circulation. Several hot‐spots of intense Fe delivery to the open ocean are maintained by standing meanders in the mean current and enhanced by transient eddies and seasonal oxygen depletion. Our results highlight the importance of fine‐scale dynamics for the transport of Fe and shelf‐derived elements from continental margins to the open ocean, and the need to improve representation of these processes in biogeochemical models used for climate studies. 

In a recent paper [Chu (2023; Chu23)], the author formulated the equations governing atmospheric motion in a spheroidal coordinate system. Since the mass distribution of the Earth is not exactly spheroidal, the true gravity is not vertical in that coordinate system. Chu23 compared the magnitude of the static horizontal component of gravity in that system to those of the dynamically active forces and concluded that the horizontal components of gravity should not be neglected. In recent papers by the authors [Chang and Wolfe (2022; CW22) and Stewart and McWilliams (2022; CW22)], we explained that the actual interpretation of the approximation made in atmospheric and oceanic modeling is not neglecting the horizontal component of the true gravity, but is a geometrical approximation, approximating nearly spheroidal geopotential surfaces with bumps on which the true gravity is vertical by exactly spheroidal surfaces. We showed that under such an interpretation, the errors due to the geometrical approximation are small. Chu23 claimed that CW22 and SM22 erroneously neglected the gravity perturbations in their analyses. Here, we explain further the differences between these approaches, in the process showing that the criticisms of Chu23 on CW22 and SM22 are invalid, further supporting our conclusion that the horizontal component of the true gravity is not relevant in ocean and atmospheric dynamics. Physically, the reason why horizontal gravity is irrelevant in the coordinate system used by Chu23 is that it is balanced by a static horizontal pressure gradient force. 

Recent high-resolution large-eddy simulations (LES) of a stable atmospheric boundary layer (SBL) with mesh sizes N=(5123,10243,20483) or mesh spacings ▵=(0.78,0.39,0.2) m are analyzed. The LES solutions are judged to be converged based on the good collapse of vertical profiles of mean winds, temperature, and low-order turbulence moments, i.e., fluxes and variances, with increasing N. The largest discrepancy is in the stably stratified region above the low-level jet. Subfilter-scale (SFS) motions are extracted from the LES with N=20483 and are compared to sonic anemometer fields from the horizontal array turbulence study (HATS) and its sequel over the ocean (OHATS). The results from the simulation and observations are compared using the dimensionless resolution ratio Λw/▵f where ▵f is the filter width and Λw is a characteristic scale of the energy-containing eddies in vertical velocity. The SFS motions from the observations and LES span the ranges 0.1<Λw/▵f<20 and are in good agreement. The small, medium, and large range of Λw/▵f correspond to Reynolds-averaged Navier–Stokes (RANS), the gray zone (a.k.a. “Terra Incognita”), and fine-resolution LES. The gray zone cuts across the peak in the energy spectrum and then flux parameterizations need to be adaptive and account for partially resolved flux but also “stochastic” flux fluctuations that represent the turbulent correlation between the fluctuating rate of strain and SFS flux tensors. LES data with mesh 20483 will be made available to the research community through the web and tools provided by the Johns Hopkins University Turbulence Database. 

Abstract Surface gravity wave effects on currents (WEC) cause the emergence of Langmuir cells (LCs) in a suite of high horizontal resolution (Δx= 30 m), realistic oceanic simulations in the open ocean of central California. During large wave events, LCs develop widely but inhomogeneously, with larger vertical velocities in a deeper mixed layer. They interact with extant submesoscale currents. A 550-m horizontal spatial filter separates the signals of LCs and of submesoscale and larger-scale currents. The LCs have a strong velocity variance with small density gradient variance, while submesoscale currents are large in both. Using coarse graining, we show that WEC induces a forward cascade of kinetic energy in the upper ocean up to at least a 5-km scale. This is due to strong positive vertical Reynolds stress (in both the Eulerian and the Stokes drift energy production terms) at all resolved scales in the WEC solutions, associated with large vertical velocities. The spatial filter elucidates the role of LCs in generating the shear production on the vertical scale of Stokes drift (10 m), while submesoscale currents affect both the horizontal and vertical energy fluxes throughout the mixed layer (50–80 m). There is a slightly weaker forward cascade associated with nonhydrostatic LCs (by 13% in average) than in the hydrostatic case, but overall the simulation differences are small. A vertical mixing schemeK-profile parameterization (KPP) partially augmented by Langmuir turbulence yields wider LCs, which can lead to lower surface velocity gradients compared to solutions using the standard KPP scheme.